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a ( erovelizéd daescriois, function is derived for a dissymmetric 
relay wit D, ©, componenvus di voth the insut and output paths, Sethods 
of evaluating tie D, C,. co,monents are discussed, The describin,, funce 
tion is applied to tae analysis of ste» and disturbance inouts, 

the descrioing function is modified to allow conputation of the 
D. C, and sinusoidal co ponents orf the error signal when a ramp input 
is apolied. vhe applicaoility of the method in various situations is 
@aseussed, 

the validity of tne describing function analysis was investigated 
by comparing calculated responses with actual system resvonse, as sinu- 
dated on the analog co.outer, 
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counsel and perceotive suggestions provided by Professor G, J, Thaler as 


mcuLtLy aavisor for this project, 


iad 








Chapter 


a 


net i: 


IV 


Vv 


T..2 20DUC SLO. 


Sige ia PLOW IOF D, C, SiGeaL LEVSLS 


THe, 


41 


seneral Considerations 
D0, ©. Levels in the Outout Path 


Soe bee la ies lane ral 


DI Sse fei C RELAY SsRVO 
JLscussion of the Describing Function 


Lae wmescripine Wiuce1on- lor a Relay wan 
no ilysteresis 


me Describing Function for a Relay with 
aysteresis 


Comparison of syunetric and Dissymmetric Relays 


Verification of Proposed Describing Functions 


GSMSRALIZED DeoCRISING FULICTION FOR A RAMP IUPUT 
General Considerations 

Describing Functior. for One=Sided Cperation 
uxoeriznentel Procedures 

Verification of the Proposed Describing Function 


The Rass Describing runction for Two-Sided 
Overation 


sxtersion of the Ramp Describing function to Other 
iypes of Non-linearities 


He COM DAL TQ@E>erOR FUTLES DT VesTicaTion 


Pave 





NI nr we) ee) 


‘O 


7 


vhanter 
BIBLIOGRAPHY 


APPanDIA J 


aPPRiIDIX IT 


om; om | Somat 
ee) WN LB ree WF ba er eae - vd 
ea 8s EE CAROL Ben Ger Ome 


title Pace 


DsarivarIO. CF vos GSNSRALIZED DESCRIBING 26 
PULICIIO: FOR & DISSYAMSTRIC RELAY 


DERLTVarION OF THS LIAIT OF SINUSOIDAL OPSRATION 81 
DULGVG OPURATION OF OLY Ors SIDE 


EV 





Fieure Nunbder 


el 


Det 


3.3 


344 


be? 


369 


LIS: OF ILLUSTRATIO.:S 





Descrivtion 
Slock diagram of a servomechanism containing 


a dissyymetric non-linearity. 


Pransfer characteristics of a dissymmetric 


relay with no hysteresis, 


Generalized describing function for a dissym.- 


metric relay with no hysteresis. 


transfer characteristics of a dissymmetric 


relay with hysteresis, 


Generalized describing function for a dissym- 


netric relay with equal hysteresis zones, 


Graphical construction for solution of equations 


Seance =O senGeu coe = iene 


Generalized describing function for a dissym- 


metric relay with unequal hysteresis zones, 


Schematic diagrams for analog computer setup 


for simulation of relay servonechanism and dissym- 


metric relay. 


Nichols Plot of linear portion of simulated 


relay servomecnanisn, 


UO 
ab) 
Q 

@ 


LA) 


10 


13 


15 


21 


23 


25 


ine 
4 a 


40H” 
7" 


ef 


4.3, +4 


LE. OF ILLUS SPA sTONS 

Jesceri tion 
m, 1Cal. Brus.) recordings ol steo ressouses 
for simulated relay servonechanisn with no 


hysteresis, 


Tyoical Orusn recordings for sinulated relay 


servonecnanis:: witn hysteresis, 


rransfer characteristics of relay operatin72 


on only one side. 


veneralizged radio cdescripim~ funcrion fora 


2 


relay actuated on only one side, 


Plots of Z vs [y used in sraphical solution 
for the critical amplitude of anout av wine 
syste i cnanges from one sided operation to 


THO Sided ooeration 


Illustration of wave shape of I (t) in non-« 


sinusoidal region of rary describing function 


plot. 


iaichols Plot of linear portions of systems 


used for verification of ramo describing func. 


G10n. 


ey: 


32 


oie 


1.2 


13 


19 








4o9Q - 4,21 


i-1 





LISY OF ILLUSTRALSTONS 


Description 


.lcnols Plot of exnerinental data obtained in 


initial investication of ram describing function. 


Averace ayolitude and frequency error vs 


yariation in hysteresis. 


ryoical “rush recordings of ramp resyoonses 
obtained during exverimental verification of 


—_ 


ram. descrioing function, 


fransfer characteristics or a dissyimetric relay, 


56 


OL 


76 





fable Sumber 


ay 
oh 


nt 


jae 


LV 


Y 


VI 


We 


eee 


List OF VASLEsS 


Valuesnof <4 lor <a dissyzinetric relay with no 


hysteresis, dy + do = 5.2 volts. 


Values of Ca for a dissy.metric relay with 


hysteresis, d, + dp = 6 volts, py + po = 10 volts, 


Values of generalized describins function for 
a dissyzmetric relay with py = oi Do = fae dy = 


do = Ae voneeee 


Hoxoperitenval results for dissyznetric relay 


without hysteresis, 


uxperimental results for dissymmetric relay 


with hysteresis, 


Data obtained in experinental determination of 


describing function for a ramp input, 


Recorded and predicted limit cycle values for 


ramp inputs, 


srrors in predicted amplitude and frequency 


resulting from small relay hysteresis errors, 


18 


ae 


27 


52 





Symbol 


G(s), Gy(s), Go(s) 


u(t) 


kk 


Lio. OF SYMBOLS AND ABBREVIATIONS 


Description 
inout cousand signal of a servonechanisn 
output position of a servomechanisn 
servo system error signal 
Laplace transforn of € 
input signal to non-linear element 
d.c, component of I(t) 
amolitude of sinusoidal component of I(t) 
frequency of sinusoidal component of I(t) 
output signal fron non=linear element 
d.c. component of O(t) 
coefficient of coswt in Fourier series 
expansion of O(t) 
coefficient of sinwt in Fourier series 
expansion of O(t) 
generalized describing function for a 
non-linear element 
non=dimensionalized ramo describing function 
phase shift associated with G, or G, 
the Laplace variable 
product of factors of the form (s + p) 
Laplace transforms of linear elements in 
a servomechanism 
the unit step function 


transient response gain of a linear element 





Sy iool 


-V, 


LIS? OF SYWgOLS AED ABBREVIATIONS 
Description 
frequency resvsonse gain of a linear element 
relay outout voltage 
relay pull-in voltaze 
relay drop-out voltage 
relay pull-in angle 
relay dron=-out angle 
notes in the preceding five symbols, n = 1 
for the positive side of the relay, 


and n = 2 for the negative side, 


dy + 42 


2X 
x 


aaa 





Pe ee 
ea 


subscript déenovine a critical value (of 224, 
or D, at whicn the systen changes fron one 


node of operation to another. 


ra 





aS ay Se 
es 1) a a 


Tero UC 


4 


in receme years, tae descrivia , function method nas been widely 
used in the anelysis of servorechanis..1s containing non-linear elenents 
(1). fhe basic describinz function techniaue consists Gf Udamearizing! 
the non-linear vortion of the systen by asSuying that: 

1. che inout to tne non-linear device is a pure sine wave of 

known azgplitude, 

2. the higner har-.onics in the output of the non-linear device 

may be neszlected, 

A representative describing function is then calculated by con. 
puting the axnplitude gain and ohase shift between tne input and the 
fundanental frequency of the fourier series exvansion of the output 
wave, <rhese techniques ere thoroughly discussed in standard texts (2). 

In sOme syste is, nowever, the presence of the non-linearity 
may cause a D. C. cofiponent to exist in either the input or output 
Math, or both. Inis is narticularly true with dissyrmetric non-~ 
linearities (3), A siiilar situation can be caused in many symmetric 
systens by introducing a ra‘o inout. In either case, the D, C. con 
ponent nay cause a change in the nagnitude and/or phase of the funda. 
iental frequency of the outout fro the nonelinear element, and it 
May result in a Significant change in the stability characteristics 
of the syste.i, 

It is therefore advantageous to define a generalized describing 


function, based on the following assu wtions: 
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\. ae & 2 nmuliwear ele rene is a oure sine wave 
Sec: MROSEIAST, & “ure wv, C, cohnonent, 

2, sane hisher harmonics in the outout of the non-linear 
ele..ent cay be neslected. In this case, the output of 
tne nonelinear ele :ent ay be reduced to a pure sine 
wave Suverumposed on a pure J. GC. cowonent, where 
botn ter:s say be evaluated fro. the fourier series 
expansion of the outnout wave, 

Paew 7eneralized describine’ function Still consists vor se 1c 
a.2plitude gain and onase shift between the sinusoidal co:ponent 
of tne input and tne funda.ental rrequency component of the output 
wave, 
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hee DUrpesoc Of CCl1S Lepore Es ue Develoo <6 = ener wires 
seribing function for a dissymetric relay and to investigate the 
characteristics of servouecuanis.:s containing such elezents, fhe 
report includes a discussion of sethods of deternining the magnitudes 
of the D, C, co:ponents; an analysis of system stability for step or 


disturbance inputs; and an analysis and discussion of the syste. 


response to step and ramp inouts, 
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2.1 Gefieral Considerations 
In all cases invesvizated in this report, it will be assumed 
that the olock diagram of the servo.;echanis. can be reduced to the 
form shorm in figure 2.1, vnere G,(s) and G(s) represent the linear 
portions of a tyoe 1 Servo cechanism, and ay is the linearized describing 


function which represents the non-linear element in the systen, 





Figure 2,1 3lock diagram of a servomechanisn containing a dissymetric 


non-linearity, 


In addition, it will ce assuued that the servo loop acts suf- 
ficiently like a low pass filter so that only the first harmonic and 
the D, C. component of O(t) are transmitted around the loop, There- 
fore, for steady state operation with sinusoidal, step, ra:zp, or dis- 
turbance inputs, we can approximate the inout to the non-linearity 
with equation 2.1; and the effective portion of the output fron the 
non-linearity can be aporoximated by either equation 2.2 or by 2qua- 


mlOnS 2.3 LNTOven 2). 


Qo 


So? 1G) Ss sw 


2,2 O(t) = “2 * 84COS Wt + D4Sin wt 
a 

Pe) es + Csin(wt + ¢) 

Coe = ane + oe 


Ls A tanzt =a 
1 
% will also be assumed that G,(s) represents a compensator 
having a transfer function of the for expressed in equation 2,6; and 
what Go(s) represents a rotor-load combination having a transfer func- 
tion of the fori exoressed in equation 2.7, 


250 U4 (s) = Ky INS 24) 


TT (s + 03) 





where 2, f O and os # 0 
Wis + Z4) 


2a CAS) So. eee 
é é sMi(s +p.) 


where 2, # 0 and Ds #0 


The basic Laplace transform equations for such a system are: 


fips s = 6, — 6, 
2.9 I(t) = G4(s)(, - 9) 


2.10 ©, = Go(s)0(t) 
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in tie case of Sih Meiscuroanes, rahip, of sinusoidal ianoucs, 
the D, ©, coioonent of Lhe output fro: tne non-linear element depends 
only on the linear gain of Go(s) and the applied input function, and 
is relatively easy to evaluate, 

The equation for the Laplace transfora O(t) can be obtained 
byeusine only the transfer functions of the linear portion of the 
System as in equation 2,17. Since tne steady state sinusoidal con. 
Bemgents have no effect on the D, ©. commonents in the linear portion of 
waewoyStvemy the principle of suverposition can be used to eliminate the 
transfor:is of these components when solving for the D, C, cozmonents 
only. The final value theore2 can then be avplied to obtain equation 

ee Lncrwmerrmen 


2.12, where O(t)y, Gas and I(t), represent the Laplace transforns of 


the non-sinusoidal components of O(t), ht and I(t), respectively, 


a rae fe y 
eae T(t 
Bate =o <= S10 Gy, = S 6. = I(t)q 
S=0 S=0 


In order to siiplify the following calculations, the limits of 
G,(s) and G(s) as S approaches zero can be partially evaluated to 
get equations 2.13 and 2.1''. These results can then be subst. a! 


Bieercquation 2yiz to get equation 2.15. 
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“or & Step iirout, or disturvence, of Nacnicude a, the Laplace 
transior of @, is siven vy equation 2.1.. For a steady state analysis, 
une initial tive is immaterial, and consequently the Laplace transforn 
of I(t), can be represented by equation 2,17 and the D. C, coaponent of 


C(t) can therefore be covputed by equation 2.1%. 


2.16 ©O, = Au(t) =4 
4 Ss 
a - 
2.17 I(t)g =+ 
mre. ke oe 
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SO 

H Similar argument can be applied to Sinusoidal inouts, since 
in this case, ©, is zero, 

For ramp inouts of magnitude 3, the Laplace transforn of the 
inout is given by equation 2.19, and the D, C, component of O(t) can 


be evaluated by equation 2,20, 
Ce 





2.19 9, = Biu(t) = 
2 
Ss 
a C r B 
2.20 2 =— S oe om a = ano 
ra KY gf Sk } 
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For an n*“ ortcor Swste), Lne s vers in the denoninator of 


C) 


Go(s) would be changed to s!*, Iquacvaon 2,21 can be derived for such 


a syste: by following a lire of reasoning sicilar to that used for the 
first order systen, JIvaluation of equation 2,21 snows that the d.c, 
component of O(t) will be zero for both ste> and ramp inouts for any 


systen of order higher than one, 
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23 D. C, Levels in the Inout Path 

In the preceding discussion, there was no limitation on the 
nature of the non-linearity renresented oy Gqs Since the D. ©, com» 
ponent in the output path is independent of G,. In ace however, 
the D, C, level in the inout oath will not only depend on the nature 
Oieune System input and the characteristics of the linear vorvion of 
the system, but will also depend on the nature of the non-linear ele- 
ment, 

in many cases the solution for the ), C, component in the inout 
Bath can be extremely cowrmlicated, and nay result in transcendental 
equations which can best be solved by graphical ~eans (3), 

it should be noted, however, that in many cases it is unnecessary 
wemcolve Gxplicitly for the 0, ©, level in the inout rath in order to 
obtain the describing function curve, This fact is particularly well 
illustrated in the analysis of the relay servo with a ras» inout, given 
in section 4.2, In other cases, such as the analyses for step inouts 
given sections 3,2 and 3,3, the solution for the D. ©. level indicates 


that for frequency resnonses analysis, the 2, C, level in the input path 
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3.1 Discussion of the Describing Function, 
sopendix I convains the derivation of the generalized describing 
miaevion Tor a dissyacetric relay witi the following chardeveriscics: 
1. she pull-in voltae for a positive inovt is not equal to 
the pullein voltage for a negative input, 
2. the drop-out voltage for a vositive input is not equal to 
the drop-out voltage for a negative input. 
3. Ine hysteresis for a positive inout is not necessarily equal 


Pome MyoucresiS Tors me slave anouUL, 


na 


The positive relay outout is not equal to the negative re- 
lay output, 

The inout vs, outout characteristics of such a relay are show 
Poeeteure [-1 of Avpendix I. 

Since most practical relays have vositive and negative outputs 
of equal magnitudes, no attenot was wade to analyze a relay with ure 
Sek Cutouts, fherefore, in tne remainder of this report, it will be 
assumed that the relay output voltages are of equal marnitude, 


3.2 Tne Describing Function for a Relay with no 





the input vs, output characteristics for a dissymmetric relay 
with no hysteresis can be represented oy figure 3.1. 

Equation 2.13 shows that for a servomechanisn of the tyne shor 
ie tigure 2.1, tne D0, C, component of the nelay outout 1s Zerouaen 
steady state operation after the apolication of a sinusoidal, Step, or 
disturbance input. Combining this result with equation I-15 of Appendix 


I gives; 





wt: 17+, 


wt: 217 -B5 


Figure 3,1 Transfer characteristics of a dissymmetric relay with no 
hysteresis, 
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The anvzle relationsnios of equations I-) and I-11 of Appendix I 
can now be introduced to solve for the 4, ©, coionent of the invut 
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Wave end to solve for B, and. ‘ee aS S201), Delo, 


Cyn W do + W 
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s BY {> 


for a relay wit no u'’stveresis the eixoression for Ay, as derived 


Gaeooerdix I, is 
Ce c.4 = 


mie expression for Opis eiven by equation I-16 of Anoendix I, 


Co@binins chis equation ttn eauation 7,2 gives 
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Table I Yalues of Ga for a relay with no hysteresis, 


dy tdy = 5.2 


squation 3,9 anolies only for cases in which the relay operates 


O11 both sides during each cycle, It is not recessary to consider one~ 
sided oneration @v tiis woint, however, Since cnis type of operation 
requires tnat the D, ©, level in the relay outout .mwst be different fron 
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relay QULD=LMMMSL be zero Tor teoinvsoial, Sees, Or idisrurbaice airs, 
if the D, vo, level in vhe outout is allowed to be different fro 
zero, aS in the case or a ram inout, then bovh one-sided and tio. 
srciec, Operation 1S @ossizjle, one subject of reso inouts is treaved in 
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can be represented by fisure 3.3. 


enievuLOon for a Relay wath “ysteresis, 


page a 6 
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racteriscvics for a cdissyimetric relay 


in tnis case, the 


combined results of equation 2,13 from Section 2,2, and equation I-12 


Prom Appendix I gives: 


S26) es: +1 a = a 
Bee 3 aah ; x 8, ey en) 
Sree A, rio ¢ = 6 + & 


Combining the anele relationshins given by equations I-5 througn 


1-11 with equation I-13 of Apoendix I, we obtain equation 3,12. 


if the hysteresis 
a. 


follows: 


= Oo anc the anzle 
2 = in 


Stee oi, Cees 


x 
is the saue for both sides of the relay, 


relationshins can be evaluated as 
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wt = X, 





wtz7-B 


wt 2 1 + XX, 


Wwt=27 - B, 


Figure 3,3 Transfer characteristics of a dissyrmetric relay with 
hysteresis, , 
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3.7 sin® - sind, = 2sini (A, = X, oos$(F, + x, ) 


a4 ee = ak 
2sin3 (8, X,,) co 38, 1 x) 


By coabiiiing the results of equations 3.11, 3.13, and 3.14, it 
can be show that X, = %, and &, =i sane D, GC, component of the 
iiyout can then be evaluated by equation 3.16, and the angles can be 


evaluated oy equations 3.17 and 3.15. 
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fhe general expression for be is given by equation I-14 of 
mopendix I, Evaluating this equation by inserting the results from 
equations 3.17 and 3.15 gives equation 3,19. The describing function 
equation, 3.20, is then obtained by substituting equations 3,12 and 
Bei? into equation 3.0, The describing function in this instance also 
reduces to exactly the sane form as the describing function for the cor- 


responding syrmetric relay, 
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rable Il contains the values of the generalized describing func- 
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EioertoOr 3 relay with dy so oe oe oie a! D4 ars 10, <A graph of this 


describing function is shoim in figure 3.4, 
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Table If Values of G, for a relay with hysteresis, 
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Ld ol [OUSTS uF oS ei sa.ie for bocli sides of the relay, 
the transcomlental eqweétiotis si.ici: must be solved in orcer to determine 
the 9, C, cotMonent of I(t) vecaie mach ::ore cormlex, They are difficult 
to solve With trizono:retric relationsnivs alone, but it is relatively 
easy to solve thea by srashnical .evnods using the construction shown in 
fisure 3,0, “his construcuion can be developed ir. the following manner: 

1, velected a nrovractor with radius &, 

2, Calculate the followins quantities for each value of X to be 


investircacec, 


pig er pea ie ae 

| ies 3.5 ane 
eg +e 

9° 

o = 22 D= 425 


3%. Construct ficure 3.5 using the calculated diniensions, 
4, ove the nrotractor along the vertical axis until 
B, + os =, ++ = as showm on figure 3.6, This satisfies 
equation 3.11. 

The construction suarantees that equations Ie3 through I-11 are 
satisfied if and only %, A, %, and 4, are selected as shown in figure 
3.6, and the D, C, cowonent of I(t) is established by the center of the 
protractor, Vicure 3,6 therefore represents the desired solution for the 
angles, and «= =, 

ane generalized describing function can then be evaluated using 
equations 3.0, 3.12, anc Iei1+, table IIT gives the values for such a 
describing function for a relay with a; or Po = 7; d, =e ane do = LL, 


this cescribing function is vlotted in figure 3.7. 





je 
2 
KS 
@ 
8 
\ 


R(P-W) 


——— 


— ms 
Rd, a W) 
3 8, R(da+W) 


| X 


. 6 


ray 

5 
mi 
® 
> 


Figures 3.5 and 3.6 Gzarhical construction for solution of equations 
3.11 and I-8 through’ I-11, 








SAI ong kD Tipe nt Ne a en ey a tO Ne i TT pa ao CN, i eRe Ah? I ep Ramat Re a’ 


| 
a 


Pe ai i TE Rl a le ee eri aa 











=1,800 . 1290 21.4 
~1,615 Move 
~1, 962 noes 10,2 16225 
~1, 4:0 rere SPS 5 ea 
| 16,0 aoe 4164 nS omy, 
: 2.0 #1,520 £0999 Sai 20.0 
= ae eG ts M60 ee 














POP, se 


Values for Gg and » for larzer valves of X are approximately 





whe savie as correspondins values in table II, | 


iadle IIL Values of the generalized relay describing function for 
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= 2 Po = 7, dl, = 2, andd, =4, 
3.4% Gonparison of Svmetric and Dissymmetric Relays, 

Equations 3.5 and 3,20 indicate that the generalized describing 
funcvtions for dissywetric rel¢ys are exactly the same as the corresvonding 
describing functions for syrmetric relays, vrovided the dissyumetric relays 
have equal hysteresis on each side, A servomechanism containing such a 
dissymuctric relay would therefore resoond to @ sinusoidal, step, or 
disturoance input in a manner quite similar to a servo containing a syne 
metric relay. 

if the servo response results in a limit cycle or sinusoidal out- 
put, the output oscillations for the syrmetrical relay would be centered 
about the position covzanded by the D, C. portion of ©., while the voutour 


oscillations for the dissy metrical relay would be offset slightly to 
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llss7 ™evtric relazy coes not have equal nysteresis zones on 
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eae7 side, une generalized dcescrioing function wey we quite daiterenc 
frox the describing functicn which vould be obvainec by assuring equal 
hysteresis zones and using equation 3.20, i. comparison of figures 3.4 
and 3,7 indicates that these differences will be relatively ninor in most 
oractical cases, and that adequate results can be obtained by using equa- 
Pods 3,20, 

Peo erilicetio. of=Fronosed Describing Functions, 

In order to verify the describing functions »reviously derived, a 
Geetoton convrol Syste? consisting of the dissymmetric relay and third 
Orcer line@r coiponents was postulated, Linear transfer functions witna 
typically low pass characteristics were selected for this systen, 

Ficure 3,0 is a scheiatic of the control system together with the 
analog corputer Simulation, this simulation provides for operation of 
the relay under varying derrees of dissyrmetry and hysteresis, (4) 
Figure 2,9 is a Nichols Plot of the linear transfer functions with an 
PEeierery gain factor of 20, 

Dee | Dissyirietric Seley <Athout ilysteresis. 

fhe relay was overated without hysteresis and with a constant 
total dead zone of 5,1 volts, vhe vlus and minus vull-in voltages were 
varied on a series of runs to obtain various conditions of dissymetry, 
step inouts of 10 volts and 15 volts were applied to the system, The 
gain for the linear portion of the system was adjusted to different 


values so as to vlace the linear curve in the desired position on the 
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fable IV is a susiary of the results obtained, Typical Grush re- 
cordings for several values of sysven eain and inout magnitudes are dis- 


played in fipures 3,10 throuyzh 3.12, which illustrate the vredicted steady 


state d.c, error levels and sinusoical characteristics, 
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eee 2 Dissyrrictric nelay with rsteresis, 
The relay was onevate! with the following characteristics: oull 
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meen) VOLtS anc - 7,35 volts; drop out, + 4.0 volts and = 1,95 volts, 
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relay Witn Krstercsis, 2 SMavid sae “ower. Lhac 2 small difference pettee: 
whe observed afd preaicted valiics is “o be emoected Since the exocrimenzval 
relay characterisvics wore noc emactly natched to those used for predic- 
tions, -.e predicted valves were based on figures 9.4 teres jie 95 nl 
and d4 + C9 = 6, che characteristics of the exmerizental relay were: 
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wen a relay servo is sWiejected to a ste» or disturoance in out, 


Y 


the servo will civucr co 2 to rest while the error remains within the 
tees Of the dead zone, or itwill enter a liseit cycle which operates 
Seewesioe of tne relay durimg eacn cvele, voimilarly, when a sine inpuc 
is applied, tne steady state error sicnal will either vary within the 
[ewes Ol une dead some, or it ill cause the relay to operate on both 
Siges durinc each cycle, 

macn Suojected to a rai indat, however, the error will never re- 
Meeeieiii tie acad Zone. she servo ouvpul Bay came to rest, but the in- 
ereasin: rano Sivcral in che inout channel imll always cause the relay 
to actuate acain, sherefore, the servo .ust enter a linlt cycle when- 
Peewee at iiout 1S anplied, imless tne ram sicnal is so large that 
tne relay does not droym out. rhis lit cycle way involve actuation 
of both sides of the rela,;, or it \1a;7 involve only tne reneated operation 
of one side, 

fre case in which tne relay overates on only one Side is the 
Simplest, and will be discussed in Section 4.2, ‘men the relay overates 
on both sides, the solwcion :wst oe obtainec. by sraphical .iethods, as 
explained in Section 4,3. 
%2 Describing Function for One-sided Operation. 

When the relay ooerates on only one side, the input vs. outout 


characteristics can be vorcrayec vy fisure “.1. ‘Yhe fourier exsansion of 


the output weve, waicn is derived in Asoondix I, becokes 
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wt= 7-6 


Figure 4,1 Transfer characteristics of relay operating on only one side, 
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2 0, = Ef (cosB + cos a) 


@e1. the servo 1.as a ram inpus, o, = Stu(t), the D. ¢, component 
Smee relays Ougsy. Bave 1S “liven by equation 2.20, “iis equacion, can 


es 


be co bined with ecauetion 4,1 to vet equation 4,4 


2.20 202 = |B 
2 Ry 


The equations for the generalized describing function can now pe 
@eraered in the following manner, without solving explicitly for the D, 


GC, component of the relay input, 
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ne MS B0inb, ih 1S converient to introduce the trigonometric 


identity siven by equation 4.6, 
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The generalized describing function for the rai input during 


one=sided operation becomes 





The relay will operate on only one side as long as X is less 


than a critical value, %,, defined by equation 4.12, 
elo ee 


Wei2 x 
; Ke 1+ sin 





if the relay has ro hysteresis, there should be no discontinuity 
in the describing function as the servo shifts from one-sided operation 
to two-sided operation, On the other hand, if the relay has a signifi- 
cant anoint of hysteresis, there will, in general, be a discontinuity 


when the shift occurs, 
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she eX feSSi0N sor vie magnitude of Gq contains no miniaua linit 
for the aiwlituce of the sinusoidal inputs The phase angle expression, 


however, indicates tuat the deseribing function is valid only if 


ee as 


P? eens om 
25 ines ll 

ZT r 

Vo 





ote that if the denominator of this expression is equal to its maxim: 
value, 2, this restriction -erely requires that the neak to peak oscil- 
lations of the Sine wave ust extend fron one side of ane hysteresis zone 
Torcuue other, 

For relays with absolutely no hysteresis, the describing function 
predicts absolutely no phase shift due to the relay. However, if the 
relay contains at least an infinitesimal a:ount of hysteresis, as all 
real relays do, there will be some small value of % for which the magni- 
tude of Ga is defined, and for which the pnase shift will be 90 degrees, 
Consequently, the effect of hysteresis must be considered to be much 
more Significant in the case of ramp inouts than for other tyves of inputs, 

The describing function defined by equation 4.11 can be vlotted 


in non-dimensionalized fora by defining the following quantities 
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fhe critical equation for stability becoines 


W412 Gaiy(s)in(s) = 3.,,5(s) 


ime plot of sla LZ! can tnen be used as the ram describing 





function for one-sided operation of any relay by simoly adjusting tne 





Saar OL the linear vlot by the factor yt == ies) oan (On ne wde= 
Do, = d 
~4 1 


scribing function of eouation “4,144 is plotted in figure 4.2, 


She aximui value of Z for which the describing function will be 





valid is 
4.19 : eee 
: amen | cre “= GUT #°sin a) 
Phe 
1 ae oC 


the describing function plot can also contain loci of equal 
values of Le Z(1 + sink). These loci will pass through points which 
correspond to the ax ai value of Z for which one-sided operation will 
occur for a relay with the corresponding ratio of dead zone to hysteresis, 
ia order to opvain tnese loci, it 1s Most convenient to first ioe 


Zvs Z(1 + sind) for various ratios of . unis requires solving 





~ Myo 
for the angle ~ using equation 4.15 then solvinz for the angle & using 
equations 4.4 and 4.7. cthe values of 2 can then be obtained for constant 


values of 2, and vhese noints car vren be plotted on the describing function 
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Figure 4,5 Illustration of wave shape of I(t) in nonesinusoidal region 
of ramp describing function plot. 
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curve, 71 Geer =, : MOS B1sés°9F 2 vs Zo + sine) for various 


> =) ° - : . ; 24 S| ° 
values oi = Se nioreetio': o8§flliese curves was wsed to obtain the 


mee 8 
=— rb 
es 
L.. 


curve for », = 10, wWiiich is shew on fikure 4.2, 
it can also oe sioim taat there is a liwiving value of phase 


shift, beyond which the relay input can no lonzer oe sinusoidal. If 


lV 


ett ,» Sinusoidal operation is possiole only if 
2 





t= 


Te 8B S .5YK. , sinusoidal operation is possible only if 


These liviting valves are also plotted on figure 4.2, The mathematical 
oroof of the lizitation is contained in Appendix IT, 


the nonesinusoidal cnaracteristics become sore and more pronounced 


wer V2 


as approaches zero or unity, and as # approaches 90 degrees, for 





example, for B= .1 Yk, and ff = 81°, the solution of equations 4,4 and 
2 

Weep requires that A = 172° and B= O°, This requires a wave shape of 
the form illustrated in ficure 4.5, The corresponding sinusoidal wave is 
srotm in dotted lines for cozmarison, 
4.3 Experinental Procedures 

As a first step toward the determination of a describing function 
for raao inputs, a series of investigative runs were made on analog 
control systems. for this purpose two control systems were simulated with 


‘a 7 Ya s = rie 
transfer functions SS and 1s » resvectively,. .ichols 


Sto + 4). ce fie a aa 


Plots for these systems are shown in figure 4,6, 
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analysis of tne results of these srelininary runs together with 


~ 


he theoretical conceots which had been developed up to this point, 





indicated the significance of 2 naraieter formed by tne ratic _B ; 
“V5 
an experinent was then designed to nore completely exolore this feature, 
B 


A series of values for the ratio , ranging from 0,01 to 0.90, 





YKy > 
were selected, Several runs were made at each of the selected values 
of the ratio, by varying B and Rn to keep the ratio constant, The 
relay was operated at a constant value of hysteresis and data was 
recorded for both the second and third order systems, The reduced data 
is shown in table VI, and on a Nichols Plot in figure 4.7, 


An examination of the data plotted in figure 4.7 leads to the 


formulation of certain conclusions, Since each set of limit cycle condi. 





tions, plotted for a constant value of ~ , seens to lie on a single 
curve it appears that a describing eee toe oor exists for each value 
of this ratio. The fact that the limit cycle conditions obtained for 
the second order system substantiate the pattem established by the 
third order system even siore clearly illustrates the existence of a 
describing function, The nature of these describing function curves ine 
dicates that they will extend across all phase angles from -90°9 to eB 
if the curves are extrapolated beyond the region of plotted data, This 
agrees with intuitive reasoning that a limit cycle must exist for all 
tyoe I systems with ramp inputs, unless the system is incapable of fole 


lowing the ramp. 
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rata obvained in experinental determination of describing 


function for ramp input. See figure 4.7 for plot, 





' + 1 
} 
fy]: 


Pal L Val Luss o f Gat 
















suaie Wy &| la ‘ 
-4 pt g Ftd Boer 4 - a Se ‘ a He eA ; 
ratahatal ae fia fe Eat letat OS) Toa Air Go am Bi as tattoo tees 
35s Gamo POBOR att eesartatats mage yo pane naam & os aati ca tiet 
: rE oa | 
i? eee A ea 
Bit) rune a ig beh 7 tr een Hi Hi ; 3 Tee ey eect yr lat 

Ree y See a | PETE Tyee i a 

| 3 : 3 | an : : | | a yeoaiteati pati FEE FA esi ct : Ht ea ‘ a t | | 

HH ey peal pes Pane 

Tap ath | ut HES 
sei be 25 ae | 
| ! H 


Loe Sag 
-. 


fees Nac peace Ba 
Baloue eee 









4 




























a 
fal 


He 
a 
ESSs 
a 
ee 

EE 


FA 
FE 
+ 
gpzsee 
His 
LEIS 
‘EEE 
S$ 
Pat 
ce 





CECE EG EE 


ees 
He 
iH snr 
ales 
Feed fieeit 
Eee 
VY 


sand sus 
aaa : 
gee 


EES 


ib 


Gate 
fase 
sans 
tH 
TH 
Aces 
His 
HE Ot 


Ee 
rus 
rH 
seater 
Bit 
fa 











nee posed Ht EA a 

ett Es EE EEREL 

ee Lat fie teen at tat 
ier 


fest 
Ee 
HH 


! ' 
a 
| 4 


sruae 
= 
uae 
sae 



















cen Haat lec eee a 

fe SEE Efe Scant eae 9 tie eta 

a Hibs Hamu 

ail iate ilies Hel rie ee 
See 


PELE Si tatestecttoatas aecapene nae 


x ee ree ‘ 5 
LI —§= t- ~ “ a i 
TI | 
4 74 dE Te \ 
| 1 
see beseeuue HA “f] : ; : ! 
oened Oo HH 4 tH Ae a es | 4 
saseifostesteattant eaansece 
1 . 
aes. ++ + i 





HE 
HEE 
: HE 
seituait 






ee 
HE 
Hea 
HEE 
HE 
i 
a 
i 
St 


SECRET ATS Sayyed es ete at ee 
2 ee eee 


sictsesezas epee ba ye Bye be 
Henle sl tl 























EEE SE Te cere ct at eee ae ter ee 

[ee satin: [FE SISIEEEISOIN gt ged (tT EEL eee HE 
scan eels eat nvr eet A 
sige ET H sregiiiit HERI HecesHeentasatititeciE ate ——o 
seater anaes 

EEHuatl riiet E eee eel ae ete et at 





Another ii: srestinz fact observed from the vlots in figure 4,7 is 


symmetry about tae value _ = 0,9), shus vhe curve for a = 0,4 
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ys td ‘vo 








coincides with that for 0.6, the curve for 0,3 coincides with that for 
0.7 etc, 
4.4 Verification of the Proposed Describins Function 

Although the exveriuient2l results presented in section 4.3 con- 
firned the existence and to some degree the nature of the describing 
function, another exoeriment was designed to explicitly verify the 
theoretically derived describing function, In addition to the second 
and third order control] systeris previously sinulated, a, fourth order 
system with transfer Function sameeren Ges used. in cone 

S(S +1)(S + 2)(S + 4) 

ducting these runs the gain of each linear system was held constant so 
that the linear curves could be superim»oosed on the plot of describing 
function curves, and the predicted linit cycles were obtained from the 
intersections, The value of the ramp inout was varied in order to obtain 


daiferent values of the ratio e the hysteresis was set at anproxie 





{y5 

mately one volt and held constant througnout, lichols Plots of the 
three linear systeris are shown in figure 4,5 and the nondimensionalized 
describing function curves are showm in figure 4,2, The predicted and 
observed limit cycles for each syster are coxupared in table VII, The 
Prush recordinss are displayed in fieures 4,9 through 4.21. 

Since some of the nredicted limit cycle values showed dis- 
erepancies of 253 or more frou those observed, it is felt that a dis- 


cussion of error sources is warranted here, First of all the inherent 


errors of any ceseribing function analysis are present, that is the 
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descrioiny Panctiom wrecicts tne ratio of the fundanental sinuscideal 
conmonents of the invut and outout simals, In many of the runs cone 
ducted (large and small ratios of = ) the input to the non-linear 
element was decidely non=Sinuscidil, ‘This can be appreciated by 
examining the S5rush recordings. It should be noted that the value 
recordec for the exoerimental value of A was simply + the peak to veak 
amplitude of the relay input, while the value of & predicted by the 
descriding function is the aimlitude of tne fundanental frequency 
component, Appendix II shows that for values of -& near 0,1 or 0.9, 
Ling 

and for values of ~ near 90°, the pullin angle approaches 1809, thus 
requiring an extremely non-sinusoidal inout signal with a peak to peak 
amplitude which approaches 3 the arolitude of the fundamental frequency 
component as illustrated oy figure 4.5, 

in section 4,2 and appendix Ii the limits for sinusoidal opera. 
tion are discussed and the area defining these limits is shown on the 
describing function curves in figure 4.2, It can be seen from the table 
of predicted and observed values tnat the errors are largest outside 
the area of sinusoidal operation, and that the predicted values become 
increasingly more accurate as the area of sinusoidal operation is ap~ 
proached, 

smother consideration in evaluating the discrepancies between pre- 
dicted and observed values is the sensitivity of the describing function 
to hysteresis, ‘The analog Sinulation of the relay employs tvio d.c, 
ayolifiers in the analog costputer, set to high gain values, It was not 


unco.mion for the drift in these amplifiers to cause changes in the 


hysteresis settince of several tenths of a volt, In order to gain an 


Ln 
are) 





appreciacion Lor “adg sarfiatu.e of errors introduced by siall changes in 


hysteresis, a revresentative system containing fourth order linear com 


“i> 


St 


nonents was exeained, The results of this cxamination are shown in 
table VIII and craphically dis»vlaved in fisure 4.8, The extrene sensi- 
tivity of the arplitude to hysteresis changes is evident, It is also 
noteworthy that these errors are essentially invariant with the ratio 
am 
ee ° 

In conclusion, after careful analysis of the data recorded, it is 
apparent that the general validity of the proposed describing function is 
proven, Specifically, when intersection with a describing function curve 
occurs within tne area defined for sinusoidal operation the amplitude 
and frequency of the limit cycle are predictable within 10; error, ‘hen 
the intersection occurs outside cf this area proportionately greater 
Gmrens are expected, but reasonable and useful estimates can be obtained, 
4.5 the Racxm Describine Function for Two-sided Qoeration, 

if the sinusoidal component of the relay input is greater than 
the critical value of X, defined by equation 4.12 or 4.19, the systen 
will enter a limit cycle in whicn both sides of the relay are actuated 
during cach cycle, In this case, the pullein and dropeout angles must 
satisfy equation 4,20, which is obtained from the sinultaneous solution 
Speeamarions 2.20 and 1-12, 

4,22 t) +8, ~of, - 8, = Sele 
Z rm 
in order to evaluate the describing function for a relay with 


hysteresis it 1S necessary to obtain the simultaneous solution of equa- 


tion 4,22 and equations I-3 throurh I-11, This solution can be obtained 
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Table VIII, srrors in vredicted ariplitude and frequency resulting from 
sall relay hysteresis errors, Zero error condition is 
for ped = 1.0, Data shown is for a system with fourth 


order linear components: 
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cra wrcally Syaaplentaie of = Py wey - FY vs, “Wdor gonstant values 





of X, Yne vaiue c 7] car teen fe sobhined fer any conbination of 
a and «, this wifomiaticn san in turn ce used to caiculate de- 
aay 


serfbins function curves wiich will be valid for twoesided oneration, 
Unfortunately, this metho: is net easily adapted to a nonedimension~= 
alized forn. Consequently it ust nernally be calculated individually 
for each relay to be investicated, 

The solution is slightiy less complicated for a relay without 
hysteresis, The transfer characteristics for such a relay are shown in 
fre sete In this case, x = B and x = ae Therefore the systen 


of Simultaneous equations which nust be soived in order to evaluate W 


consists of only equation 4,23 and equations I=9 and I~11. 


Lae on. 3 (i 


The equations for the coefficients of the fundamental components 


of the fourier series expansion of O(t) reduce to 


ht a ¢ = 0 
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The solution can be simplifisu further sy noting that 


426 tos By cos Bey = 2 cos3(8, + 8 2)cosi( F - B,) 
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427 oF = zec0s 3( B. = 8 .)cos BT 








an Gver siveler sol *bavivwen be omeined for a system containing 
an ideal relay, in this case tnere is no dead zone, and A, = Spe 
Therefore equation 4.27 reduces to eauation 4,28 and the describing 


function can be expressed by equation 4.29, 
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It should be noted that there will se a discontinuity between the 
describing function for oneesided operation and the describing function 
for two-sided overation wnenever there is a Significant hysteresis effect, 
Therefore, if tne system is operating on only one side, and a small 
change in inout or gain causes the error signal to change so that 
x a4 OF the syster outout may change radically. 

it is possible that this discontinuity may help to predict sone 
types of "jumo resonance” found in relay servos, Suppose, for instance, 
that the system is subjected to a sinusoidal innvut, o. = Asin wii, such 
that Awot = 2 fRyy if Wis mucn less than the naturel frequencies 
of the system, the servo will behave as if the input signal were a 
slowly changing rawp which has a -aximen slope of AW, at Wot = O and 
decreases to zero at bit = W 

.OW Suppose the SvStem characteristics are such that the relay 
will operate on only one Side when Wt = 0, Then the initial operating 
point will be detersined bv the rasm describing function for one-sided 
operation, The operating point will tend to travel along the linear 
curve, with the instantaneous overating point being deternined by setting 


Bes AW,COS UAL. 





° 





If the linear curve crosses the critical curve of D, for the 
relay used, the oosracins »oint will eventually reach a point where one~ 
Sided operation can no longer occur, and a two-sided limit cycle will 
result. During the next nortion of the cycle, the servo will operate 
in accordance with the describing function for two-sided operation. The 


servo will return to one«Siced operation when A decreases so that 





2&2 : ‘ a (1 + singy,), where both X and ox, are determined by the 
n 


two-sided describing function, 

In general, the magnitude of the systen cutout will "jum" to 
a larger value when the relay shifts fron one-sided operation to twoe 
Sided operation, and will "junp" to a smaller value whenever the relay 
returns to one#Sided operation, due to the discontinuity between the two 
deseribing functions, There may also be a change in the frequency of 
the limit cycle when the shift occurs, Both the frequency change and 
magnitude change can easily be vredicted once the describing function 
Curves have been plotted, 
4,6 Extension of the Ran 


o Describing Function to Other tyoes of I 





linearities, 

The rarro deScribing function can be extended to apply to other 
types of nonelinearities. In general, the fourier series exoansion of 
the output wave will sive several equations of the following tynes: 

1. A relationshio for a, similar to equation I-3, 

2. Relationships for a4 and by Similar to equations I-5 and 1-7, 

3. Relationships for all angles involved, similar to equations 


I-8 throuch In44, 


ve 





sunce ewGdcicil 2) 20 Meeblies Tor any nonelinearity, it can de 
cosibineu with the fourier equivalent of aj to get a relation between 
the rao magnitude and the angles involved, similar to equation 4.4, 
This in turn can be cornbined with the angle relationships to solve for 

he d.c. component, , Once Wis evaluated, it is a relatively sinole 
matter to evaluate the equations for a4 and by and substitute these 
values into equation 3.3 to obtain the describing function. 

It should be noted, however, that the transcendental equations 
involved will become more complex as the nature of the non-linearity 
becomes more conplex, and the solution of these equations may become exe 
trenely aifficult, in addition, other types of nonelinearities may have 
different modes of operation corresponding to the one=sided and ttuo~ 
Sided modes of operation for the relay servomechanism, and there may be 
discontinuities between the describing functions for these different 
modes, 

Thus a complete analysis will require the evaluation of a separate 
group of describing function curves for each mode of operation, along 
with a set of limit curves to define the points at which the system will 


shift from one mode of operation to another, 
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During che course or this investication, several interesting 


problems were nosed which eoneared to be worthy of further study, hese 


Suggested subjects ares 


1. 


» comprehensive study of the effect of different types of 
dissyrmetry on systens containing different tices of none 
linearities, 

. Tiore exhaustive analysis of the rain describing runction 
for a relay servovechanisn im Vinien unre relay operates on 
both sides, nis analysis should include an investigation 
of the discontinuity between the describing function for one-~ 
Sided operation and the describing function for two-sided 
operation to deter.ine tmether or not these cesceribing 
funetions can be used to predict certain types of "juno 
resonance" effects, as suggested in section 4,5, 

a thorough study of the extension of the ramp describing 
function to non-linear devices other than the relay would also 
be valuable to determine the effect of rarmo inputs on systenis 
containing Such elements. An outline for the procedure for 


such an investization is given in section 4,6, 
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DERIVATION OF Tia 


& GSSERALIZEN DESCRIBING FUNCTION 
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DISSY3 METRIC RSLAY 


A, Input vs, output characteristics, 
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D, special Cases for 14 = to 


if i, = Y5 = Y, equations I-35, T=m5, and I-7 reduce to 


A 
KK, 


T-i2 a. ==(6, + 4, = @, = &j) 

= ee : ; 

te13 aq, = =(sin B e Si x sin 6, =» Sin t .) 
Tait by = a(cos B 4 + cos X 4 + cos A, + cos X 5) 


if %,; = Y5 = Y and the relay has no hysteresis, then x Fe he 


and oS Pos and equations I-12, Ie13 and I-14 reduce to 


I-15 a, = (8, - A4) 
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ay ty = ty = i and the relay operates only on the positive side, 
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writing, so that equations 1812, T#13, and I-14 reduce to 
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